Collapse of thermal activation in moderately damped Josephson junctions 
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We study switching current statistics in different moderately damped Josephson junctions: a para- 
doxical collapse of the thermal activation with increasing temperature is reported and explained by 
interplay of two conflicting consequences of thermal fluctuations, which can both assist in premature 
escape and help in retrapping back into the stationary state. We analyze the influence of dissipa- 
tion on the thermal escape by tuning the damping parameter with a gate voltage, magnetic field, 
temperature and an in-situ capacitor. 
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PACS numbers: 74.40.+k, 74.50.+r, 74.45.+C, 74.72.Hs 

Decay of mctastable states is an important process in 
many scientific areas Dissipation plays a crucial role 
in the decay dynamics. The influence of dissipation on 
thermal and quantum escape from the superconducting 
(S) to the resistive (R) state in Josephson junctions has 
been intensively studied both theoretically Q, 0, 01 
and experimentally El 0, 0, El El El j most recently in 
connection with the problem of decoherence in quantum 
systems [Tl| . 

So far switching statistics was studied for 
superconductor- insulator- superconductor (SIS) 
junctions, while superconductor- normal metal- super- 
conductor (SNS) junctions, which arc characterized 
by stronger dissipation effects, remain, to our knowl- 
edge, unstudied. Analysis of dissipation effects in SIS 
junctions is complicated by ill-defined quality factor 
which can not be represented by a simple constant 
Conflicting reports exists about what kind of resistance 
R determines the effective damping of SIS junctions: the 
normal resistance 0], the high frequency impedance of 
circuitry 0, or the quasiparticle resistance Q. This is 
not a problem for SNS junctions, which typically have a 
resistance R considerably smaller than the open space 
impedance ~ 3770, and are well described by the RCSJ 
model with frequency independent Q. 

Here we study switching current statistics in moder- 
ately damped superconductor-two dimensional electron 
gas-superconductor (S-2DEG-S) and Bi 2 Sr 2 CaCu208+5 
(Bi-2212) high-T c intrinsic Josephson junctions (IJJ's). 
Being able to tune the damping parameter by a gate volt- 
age, magnetic field, temperature and an in-situ shunting 
capacitor, we analyze the influence of dissipation on the 
thermal activation (TA). For both of the drastically dif- 
ferent systems we observe a sudden collapse of TA with 
increasing T and explain this paradoxical phenomenon by 
interplay of two conflicting consequences of thermal fluc- 
tuations, which on one hand assist in premature switch- 
ing and on the other hand help in retrapping back to the 
S-state. We present numeric and analytic calculations 



which are in good agreement with our experimental data. 

Fig. 1 shows typical Current- Voltage characteristics 
(IVC's) for a) a S-2DEG-S (Nb-InAs-Nb) junction #1 
and c) a Bi-2212 mesa structure containing nine IJJ's. It 
is seen that the IVC's are well described by the RCSJ 
model with frequency independent Q E2- Details of 
sample fabrication and characterization can be found in 
Refs. [H El and [H for S-2DEG-S and IJJ's, re- 
spectively. Measurements were done in a shielded room 
in a dilution refrigerator for S-2DEG-S, or in a He-4 cryo- 
stat for Bi-2212. Switching currents were measured using 
a sample-and-hold technique. 

From Fig. 1 it is seen that the IVC's exhibit hysteresis. 
Figs. 1 b) and d) show T— dependencies of the switch- 
ing, Is, and retrapping, Jr, currents for other Nb-InAs- 
Nb and Bi-2212 junctions studied here. According to the 
RCSJ model, Q can be obtained from the magnitude of 
hysteresis, Is/Ir- In case of Bi-2212, the experimental 
Is/Ir agrees well with the calculated Q using capaci- 
tance of IJJ's C ~ 68.5 fF/fim 2 0. For the unshunted 
S-2DEG-S #1 from Fig. 1 a) the I s /Ir ^ 1.3, would 
correspond to Q ~ 1.4 and C = 0.1 lpF, consistent with 
the estimated value of the stray capacitance. 

On the other hand, the hysteresis in SNS junctions 
can be also caused by self-heating El > non-equilibrium 
effects or frequency dependent Q 0. In order to 
understand the origin of hysteresis, we fabricated an in- 
situ shunt capacitor, consisting of AI2O3/AI double layer 
deposited right on top of the Nb-InAs-Nb junction. The 
IVC's of the S-2DEG-S #1 before and after C-shunting 
are shown in Fig. la). It is seen that the hystere- 
sis increased considerably, while R was little affected by 
C— shunting. Such behavior is inconsistent with the self- 
heating scenario. Thus the hysteresis in our junctions 
is predominantly caused by the finite Q > 1. A similar 
conclusion was made for other planar SNS junctions [l9| . 
where it was observed that there is no correlation between 
the hysteresis and the dissipation power at / = Jr. 
S-2DEG-S junctions provide a unique opportunity to 
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FIG. 1: a) IVC's of a Nb-InAs-Nb junction #1 at T = 30mK 
before and after deposition of an in-situ shunting capacitor; 
b) Measured switching (solid lines) and retrapping (dashed 
lines) currents of another Nb-InAs-Nb #3 at four magnetic 
fields, c) IVC's of a Bi-2212 mesa containing nine stacked 
IJJ's at four different T < T c ~ 93K. Multibranch struc- 
ture is due to one- by-one switching of IJJ's. d) Measured 
switching Is and retrapping Ir currents of another mesa. 



tune the Josephson coupling energy Ejq and Q by ap- 



uplmg ei 

plying gate voltage V g [13, [l4| . For this purpose a thin 
gate electrode was deposited on top of the InAs. Fig. 
3 a) shows switching current histograms at T = 37mK 
for S-2DEG-S #2 at different V g . The inset shows the 
width at the half-height of the histograms AI vs. the 
most probable switching current Ismax- It is seen that 
initially histograms are getting wider with increasing neg- 
ative V g , consistent with the increase of TA with decreas- 
ing Ejq/T. However, at V g < —0.35V a sudden change 
occurs and AI starts to rapidly collapse. 

Fig. 2 b) shows AI vs. T for the S-2DEG-S junction 
#3 at four magnetic fields, the same as in Fig. 1 b). In 
all cases we can distinguish three T— regions: 

(i) At low T the histograms are independent of 
T, indicative for the Macroscopic Quantum Tunnelling 
(MQT) regime 0, S H @. The decrease of AI with H 
leaves no doubts that we observe for the first time the 
MQT in SNS junctions (details will be published else- 
where) . 

(ii) At intermediate T, AI increases in agreement with 
the TA calculations, shown by dashed lines in Fig. 2 c), 
for which the escape rate from S to R state is given by 



Tta =o,t— exp 
Ztt 



AU 

k B T 



(1) 



Here AU is the barrier height, and io v = uj p q(1 
(I/Ico) 2 ) 1 / 4 is the Josephson plasma frequency, lu p q 



a 
s 

o 

U 




3.0 3.5 4.0 

S-2DEG-S #3 



4.5 5.0 5.5 

/(uA) 



< 
c 



unshunted' '.J.^ 1 - 



b) y. 



140 
120 
100 

60 ^ 

40 
20 




Simulation: 
C=0.15pF 
R =10 £2 



c) 



y 

■V. 

_i i 



200 400 600 800 

r(mK) 




200 400 600 800 

r(mK) 



FIG. 2: Switching statistics of S-2DEG-S junctions: a) his- 
tograms at different gate voltages for junction #2. Inset 
shows the width vs. the most probable switching current 
Ismax- A sudden collapse of A7 occurs at V g < —0.35V. b) 
The width of histograms vs. T for the junction #3 at four 
magnetic fields, the same as in Fig. lb). Three T— regions 
can be distinguished: the MQT region at low T, the TA re- 
gion at intermediate T, and collapse of histograms at higher 
T. c) Numerical simulations for the case of Fig. 2b). Dashed 
and solid lines represent AI for classical thermal activation 
disregarding and taking into account retrapping, respectively. 



(2eI c o/hC) 1 ^ 2 and I c o is the fluctuation- free critical cur- 
rent. Dissipation enters only into the prefactor of Eq.JI}, 
which for our moderately damped junctions is a t = 
(1 + 1/4Q 2 ) 1 / 2 - 1/2Q 3j, where Q = uj p RC. 

(iii) At higher T, the histograms start to rapidly col- 
lapse leading to a downturn of AI. This paradoxical 
phenomenon is the cental observation of this work. 

Fig. 3 c) shows that a similar collapse occurs in IJJ's 
at T* ~ 75K. Here we show the effective escape temper- 
ature T esc , which indicates how much the relative width 
AI/I c0 differs from the TA prediction, Eq.JTJ). Figs. 3 a) 
and b) show switching current histograms of a single IJJ 
just before and after the collapse. As reported previously 
, at T < T* the histograms are perfectly described by 
TA, shown by the dashed lines. However, at T > T* the 
histograms become narrower and loose the characteristic 
asymmetric shape, as seen from comparison with the TA 
simulation in Fig. 3b). 

We start discussing the observed phenomenon by ex- 
cluding scenarios which can not explain it. First, it can 
not be due to T— dependence of the damping parameter 
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FIG. 3: Switching histograms of a single IJJ at a) T < T* and 
b) T > T*: circles represent experimental data, dashed lines 
- TA simulations, Eq.Q, and dashed-dotted lines - proba- 
bilities of not being retrapped, EqlJ^J. Solid lines show the 
conditional probability of switching without being retrapped 
(the product of dashed and dashed-dotted lines). Note that 
both the width and the shape of the histograms change at 
T > T* c) The effective escape temperature for the same IJJ. 
A sudden collapse of T esc at T* ~ 75K is seen. 



since Q(T) changes only gradually through T* and since 
we did take into account the Q— dependence of the TA 
prefactor at in our simulations. Second, it can not be 
caused by frequency dependent damping due to shunting 
by circuitry impedance. Indeed, we also observed a simi- 
lar collapse for planar SNS junctions [19j with R ~ 0.217, 
for which such shunting plays no role. 

To explain the phenomenon we first note that T* is 
close to the temperature at which hysteresis in IVC's 
vanishes (cf. Figs. 3c, Id and 2b, lb), indicating that 
retrapping plays a role in the observed phenomenon. The 
rate of TA retrapping from R to S state is known only 
for strongly underdamped junctions Q 1 
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, (2) 



where Qo = u p oRC . Note that unlike the TA escape, the 
TA retrapping depends strongly on dissipation[lfJ, due 
to the Qq factor under the exponent in Eq.|J2J. 

The probability to measure the switching current / is 
a conditional probability of switching Ps(I) from the S 
to the R-state and not being retrapped back during the 
time of the experiment: 



PnR = 1 



C ° P R (I)dI/ [ C " P R (I)dI, (3) 



l-tf°P R (I)dI 



is the retrap- 



where P R (I) = 

ping probability. 

Dashed-dotted lines in Fig. 3 a,b) show the calculated 
P nR . It is seen that at T < T* the P nR — 1 in the 
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FIG. 4: Numerical simulations for S-2DEG-S #3: a) In- 
dependence of the fluctuation free I c o — 20fj,A (dotted line), 
the most probable switching current Ismax (solid line) and 
the most probable retrapping current lR ma x (dashed line), 
b) The width of switching histograms disregarding retrap- 
ping (dashed line) and taking into account retrapping (solid 
line), c) T esc vs. T for S-2DEG-S #2* before and after in- 
situ C-shunting. d) The height of escape barrier (at the most 
probable switching current Ismax) vs T*: symbols represent 
experimental data from Fig. 2b), the solid line corresponds 
to the simulation in Fig. 4 b). 



region where P5 > 0, therefore retrapping is insignificant. 
However, at T > T*, retrapping becomes significant at 
small currents. The resulting conditional probability of 
measuring the switching current, P(I) = Ps(I)P nR {I), 
normalized by the total number of switching events, is 
shown by the solid line in Fig. 3 b). It is seen that it 
explains very well both the reduced width and the almost 
symmetric shape of the measured histogram. 

Figs. 4 a,b) show results of simulations, in which we 
intentionally disregarded the T— dependence of I c q, keep- 
ing Qo an d Ejq constant. In the simulations we used the 
parameters of S-2DEG-S #3. It is seen that the most 
probable retrapping current I Rma x has a weak T— de- 
pendence, consistent with the experiment, see Fig. 1 
b). On the other hand, Ismax decreases approximately 
linearly with T and eventually crosses I Rma x- Fig- 4 b) 
represents the width of histograms. It is seen that within 
the TA model AJ continuously increases with T. How- 
ever, retrapping reduces AI as soon the switching and 
retrapping histograms start to overlap. 

The T* can be estimated from the system of equations: 



max 

) ~ (di/dt)/i c0 , 

^ R {Tdown) = TtA- 



(4) 
(5) 



Here the first equation is the condition for Ismax, which 
states that the junction should switch into the R-state 
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during the time of experiment. From Eas. (|ll4|) it follows 
that MJ(I Smax )/k B T ~ In 



atLJpI c 



[2^lJTt\ = F ' which a S reeS 
with experiment, as shown in Fig. 4d). Taking AU ~ 

(4v / 2/3)£'jo [1 - I S /Ico} 3/ '\ and neglecting T— depen- 
dence of I c q, we reproduce the linear T— dependence 
Ismax/Ico - 1 - (3Yk B T)/(4y/2Ejo), seen in Fig.4 a). 
Substituting this expression into Eas. 15121) and assuming 
I no — (valid for Q > 2), we obtain: 
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From Eq.© it follows that T*/Ej depends almost 
solely on Q . Fig. 4 c) shows T esc vs T for a S-2DEG- 
S #2* (similar to #2) before and after C— shunting. A 
dramatic difference in the behavior of TA is obvious. As 
shown in Fig. la) the C— shunting affects almost solely 
Qo. Therefore, switching from S to R state is not strongly 
affected by C— shunting. On the contrary, retrapping is 
affected considerably because Jrq ~ l/Qo- Under these 
circumstances, higher T is required to reduce Ismax to 
the level of Ir, resulting in the increase of T*. 

To get an insight into the phase dynamics at T > T*, 
we show in Fig. 4 d) the dependence of AU(I = Ismax) 
vs. T* for the case of Fig. 2b). It is seen that AU 
scales with T. The solid line in Fig. 4d) corresponds 
to AU(Is m ax)/k B T = 24.3 ~ Y obtained from simula- 
tions presented in Figs. 4 a,b) and demonstrates excellent 
agreement with experiment. The large value of AU /k B T 
implies that the junction can escape from S to R state 
only few times during the time of experiment. There- 
fore, the collapse is not due to transition into the phase- 
diffusion state, which may also lead to reduction of AI, 
previously observed in SIS junctions jill2 1| (even though 
we assume that the decrease of AI reported recently in 
Ref.|2H may also be described by our model). Indeed, 
phase diffusion requires repeated escape and retrapping, 
which is only possible for AU/k B T ~ 1 0, 113 Careful 
measurements of S-branches in the IVC's at T > T* did 
not reveal any signature of dc-voltage down to ~ lOnV 
for S-2DEG-S and - lfiV for IJJ's. Furthermore, as 
seen from comparison of Figs, lb, 2b) and Figs. Id, 3c), 
the IVC's remain hysteretic at T well above T*, which 
is incompatible with the phase diffusion according to the 
RCSJ model Q. As can be seen from Fig.l c) the first in- 
dication for the phase diffusion in our IJJ's appears only 
at T > 90K, meaning that all the collapse of TA shown 
in Fig. 3 c) at 75K < T < 85K occurs before entering 
into the phase diffusion state. 

For a quantitative comparison with experiment we per- 
formed full numerical simulations of Eqs. (1-3) taking 
into account the T— dependence of I c q, shown in Fig. 1 
b) and the exact value of hysteresis I c o/Iro within the 



RCSJ model. Results of the simulations for the S-2DEG- 
S #3 at four magnetic fields, corresponding to Figs, lb) 
and 2b) are shown in Fig. 2c). Dashed and solid lines 
represent the simulated width of histograms disregard- 
ing and taking into account retrapping, respectively. It 
is seen that simulations quantitatively reproduce T* for 
all four magnetic fields. The capacitance C = 0.15pF, 
which was the only fitting parameter, is the same for 
all four curves and corresponds to the expected value of 
stray capacitance. Taking into account that T* is very 
sensitive to 7 c o and C, see Eq.©, we may say that the 
agreement between theory and experiment is excellent. 

In conclusion, we observed a paradoxical collapse of 
thermal activation with increasing T in two very differ- 
ent types of Josephson junctions with moderate damping. 
The phenomenon was explained by the interplay of two 
conflicting consequences of thermal fluctuations, which 
can both assist in premature switching and help in re- 
trapping back into the S-state. The retrapping process 
is significant at small currents, causing cutting-off the 
thermal activation at small bias. We have analyzed the 
influence of dissipation on the thermal activation by tun- 
ing the damping parameter with the gate voltage, mag- 
netic field, temperature and in-situ capacitive shunting. 
Numerical simulations are in good agreement with exper- 
imental data and explain both the paradoxical collapse 
and the unusual shape of switching histograms. 
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